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ABSTRACT 

We present a numerical study of rapid, so called type III migration for Jupiter-sized planets 
embedded in a protoplanetary disc. We limit ourselves to the case of outward migration, and 
study in detail its evolution and physics, concentrating on the structure of the co-rotation and 
circumplanetary regions, and processes for stopping migration. We also consider the depen- 
dence of the migration behaviour on several key parameters. We perform this study using 
global, two-dimensional hydrodynamical simulations with adaptive mesh refinement. We find 
that the outward directed type III migration can be started if the initial conditions support 
Z > 1, that corresponds to initial value > 1.5. Unlike the inward directed migration, 
in the outward migration the migration rate increases due to the growing of the volume of 
the co-orbital region. We find the migration to be strongly dependent on the rate of the mass 
accumulation in the circumplanetary disc, leading to two possible regimes of migration, fast 
and slow. The structure of the co-orbital region and the stopping mechanism differ between 
these two regimes. 
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1 INTRODUCTION 



The discovery of extra-solar planetary sy stems I May or & Ouelozl 
ll995tlMarcv et al.ll2000l.lVogt et al.ll2002n has invoked a strong in- 
terest in the theories of planetary system formation. The occurrence 
of so-called 'hot Jupiters' (objects with minimum masses M > 
M-jj. and semi-major axes a < 0.1 AU) has shown that planets can 
be found on orbits very different from the formation sites of plan- 
etary cores. Since the in situ form ation of these objec ts is difficult 
both in the core accretion scen ario iPollack et aljl 996) and through 
direct gravitational instability (Boss 200 lh . the global migration of 
planetary cores has become the subject of a large number of stud- 
ies. The inward migration of the planets d ue to planet-disc grav- 
itational inte raction was first described bv lGoldreich & Tremaind 
dl979lll98Ch . prior to the discovery of extra-solar planetary sys- 
tems. Later on it was worked o ut in a large numbe r of inv e stigations 
summ arise d in the reviews bvllin & P apaloizou dl993h . iLin et al.l 
( 2000) and Papaloizou et alj d2007l) . 



In the standard picture the planet exchanges angular momen- 
tum with the disc through the excitation of spiral density waves at 
the Lindblad resonances and through the interaction at the corota- 
tion resonance. However until recently, the second mechanism was 
assumed to be inefficient and was therefore neglected. The non- 
axisymmetric pattern of the waves excited at Lindblad resonances 
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acts back the planet providing the torque driving its orbital migra- 
tion. This leads to two possible regimes of the migration: type I 
for planets embedded in the disc and type II for planets opening 
the gap. In both cases the direction of migration usually is in- 
ward. Note that type II is however capable of moving the planet 
outward provided the disk has local density gradien t s supporting 
outward viscous spreadin g dLin & Papaloizoj|l993l : I Ward 120031 : 
ICrida & MorbideUill 20071) . 

Recently the corotation resonance have been found to be 
important and capable of slowing down an d even reverting the 
direction of migration in typ e I regime jMasset et alj 1 20061 , 
Paardekooper & Mell emal2006t) . However, it can also lead to a new 
and very fast migration mode (called type III) that is not truly res- 
onant, but depends strongly on the gas flow in the planet vicinity 
and does not have a predetermined direction jMasset & Papaloizoul 
2003; I Artymowic z 2004, Pa paloizou et alj|2007h. This type of mi- 
gration was studied numerically by [Masset & Papal oizou d2003t) 
who performed global two-dimensional simulations of an in- 
ward and outward directe d migration. T h e inw a rd directed 
migrat ion was studied by iD'Angelo et alj J2005h ■ IPapaloizoul 
(2005) considered local shearing box simulations in 2-D, and 
lArtvmowicz & de Val Borrol <2008h in 3-D. 

This paper is the third in a series devoted to a numerical 
investigation o f type I II migration for the high-mass planets. In 
IPeplinski et alj |2008a) (henceforth Paper I) we showed how the 
applied disc model critically influences the outcome of the simula- 
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tions. We found two physically motivated corrections to the stan- 
dard disc model, both necessary to remove non-physical effects 
from the simulations. The first one is a correction of the gas accel- 
eration (due to self-gravity) which ensures that the circumplanetary 
disc moves together with the planet. The second is a modification 
of local-isothermal approximation, which allows an increase of the 
temperature inside the planet's Roche sphere and thus decreases 
the amount of gas accumulated in the planet's vicinity. Paper I con- 
tains the detailed description of the disc model, the code and the 
convergence tests. 

This method was applied to t he stud y of inward directed 
type III migration in Peplinski et al. (2008b) (henceforth Paper II). 
There we focussed on the physical aspects of type III migration 
and described in detail the mechanisms driving rapid migration. 
We also discussed the mechanisms for stopping type III migration 
(or rather transforming it to the very much slower type II regime). 

In the present paper we focus on outward directed type III mi- 
gration. Although the driving mechanism is the same for both the 
inward and outward directed migration, the existence of the dif- 
ferential Li ndblad torque breaks the symme try between both cases 
as shown in Artvmowicz & Peplin ski (2008). The analytical theory 



m lArtvmowicz & Peplinskil d2008l) also shows that unlike inward 
migration, outward migration only is allowed for a co-orbital mass 
deficit A/a larger than one. Paper II showed that inward migration 
will always slow down due to the shrinking of the size of the co- 
orbital region. For outward migration one would expect the reverse, 
i.e. a speed up of migration, allowing the planet to travel long dis- 
tances in the disc. 

The layout of the paper is as follows. Section [2] gives the nu- 
merical setup of the disc model. In Sections|3]and|4]we describe the 
results of simulations of outward migrating protoplanets and the de- 
pendency of the migration on the various parameters. In Sections|5] 
and [6] we analyse the stopping of type III migration and eccentric- 
ity evolution. Finally in Section|7]we discuss their implications for 
extrasolar planetary configurations. 



2 OVERVIEW OF THE NUMERICAL SETUP 

The full description of the adopted disc model, the numerical 
method and the boundary conditions was given in Paper I. In this 
chapter we will give only a short description of the adopted units 
and initial conditions. 

For our simulations we use the Flash hydro code version 2.3 
written by the FLASH Code Group from the Centre for Astrophys- 
ical Thermonuclear Flashes at the University of ChicagtQin 1997. 
Flash is a modular, adaptive-mesh, parallel simulation code capa- 
ble of handling general compressible flow problems. It is designed 
to allow users to configure initial and boundary conditions, change 
algorithms, and add new physics modules. We use the code in the 
pure hydrodynamical mode in two dimensions, and the adaptive 
mesh is used to achieve high resolution around the planet (4 lev- 
els of refinement). The simulations are performed on the Cartesian 
grid in the inertial reference frame. 

In the simulations we adopt non-dimensional units, where the 
sum of star and planet mass Ms + Mp represents the unit of mass. 
The time unit and the length unit are chosen to make the gravita- 
tional constant G = 1. This makes the orbital period of Keplerian 
rotation at a radius a = 1 around a unit mass body equal to 2n. 
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However, when it is necessary to convert quantities into physical 
units, we use a Solar-mass protostar Ms = Mq, a Jupiter-mass 
protoplanet Mp = Mq^_, and a length unit of 5.2AU. This makes 
the time unit equal to 11.8/27T years. 

In all simulations the grid extends from —4.0 to 4.0 in both 
directions around the star and planet mass centre. This corresponds 
to a disc region with a physical size of 20.8 AU. 

The initial surface density E profile is given by a modified 
power law: 



Einit = i/'(r c )Eo(r- c /r-o) c 



(1) 



where r c = \r — rc\ is the distance to the mass centre of the 
planet-star system, ro is a unit distance, and ij) is a function that 
allows introducing a sharp edges in the disc. To start the outward 
directed migration we introduce a density jump at the initial planet 
position (see Paper I Fig. 2). The disc mass is characterised by the 
disc to the primary mass ratio 



A*d = 



Einit(ro)7rro _ E 7 



(2) 



Ms Ms 
We investigate different density profiles by changing as from —1.5 
to 0.0, and /^d from 0.005 to 0.01. For the Minimum Mass Solar 
Nebula (MMSN) /x D = 0.00144 and a E = -3/2. 

To enforce rapid migration the planet is introduced instanta- 
neously on a circular orbit of semi-major axis equal 0.8. We per- 
formed the simulations with M P /M S equal 0.0007, 0.001 (i.e. one 
Jupiter mass, M^, for a one-solar-mass star) and 0.0013. 

The aspect ratio for the disc with respect to the star is fixed at 
h s = 0.05, whereas the circumplanetary disc aspect ratio h p can 
range form 0.4 to 0.6. 

The smoothing length of the stellar potential, r so ft, is taken to 
be 0.5. For the planet this parameter was chosen be a fraction of 
the planet's Hill radius Rh — a[Mp/(3Ms)] , and was set to 
fsoft = 0.33i?H- The corresponding size of the envelope r cnv was 
set to 0.5i? H . 

Flash has the ability to track multiple fluids. We use this fea- 
ture to investigate the evolution of the fluid in corotation region. In 
our simulations we introduce a tracer fluid that has a value of 1 in 
the (initial) corotation region (ainit — 2Rh) < r < (ainit + 2Rn), 
and zero outside it. This allows us to distinguish between the fluid 
captured by the planet in the horseshoe region and the fluid flowing 
through the corotation region. 



3 OUTWARD MIGRATION - STANDARD CASES 

Although the driving mechanism for inward and outward type III 
migration is similar, we find that the initial conditions for each 
have to differ significantly. As found in Paper II, a relatively low 
mass disc (about twice the MMSN) is sufficient to start the in- 
ward migration of a Jupiter mass planet with an initial semi-major 
axis a = 3. To start outward migration for such a planet initially 
placed at a = 0.8, we need a three times more massive disc (with a 
sharp inner edge). Another difference is that for outward migration 
the planet accumulates a much larger amount of gas as it moves 
through the disk. This rapid increase of the effective planet mass 
Mp significantly modifies the planet's orbital evolution. However, 
the evolution of Mp depends strongly on the chosen parameters 
such as hp, the scale height of the circumplanetary disc. Thicker 
(warmer) circumplanetary discs accumulate less material, and thus 
results in a different orbital evolution. They also represent two dif- 
ferent mechanisms of stopping type III migration. 
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Figure 1. Orbital evolution of the outward migrating planet. Upper left and upper right panels show the planet's semi-major axis a and eccentricity e. The 
migration rate a and the non-dimensional migration rate Z are presented in middle left and middle right panels. The bottom row shows the torque exerted by 
the gas on the planet T as a function of time (left plot) and the non-dimensional migration rate (right plot). Curves 1 and 2 correspond to the models with a 
time-dependent Mp* = Mp = Mp + M ao f t (model Ml) and constant Mp* = Mp = 0.001 planet's mass (model M2) respectively. To remove strong 
short-period oscillations the data of Ml and M2 were averaged over 5 and 10 orbits respectively. 



To be able to compare the cases of high and low mass accu- 
mulation, we will here keep h p constant at 0.4, but instead in one 
case allow the effective planet mass to grow, by adding the mass 
content within the smoothing length (Mp = Mp, model Ml) and 
in the other keep the effective planet mass constant at the initial 
value (Mp = M P = 0.001, model M2). For both models the ini- 
tial surface density profile has fiu = 0.01 and as = —0.5, with 
an inner edge placed at r = 0.8. At this edge the surface density 
drops by a factor of 10 (see Paper I, Fig 2). 



Note that these models were used in Paper I, Sect. 5.2.2 in 
the discussion of the dependency of the planet's orbital evolution 
on the mass accumulation in the circumplanetary disc. Here we 
discuss them in more detail^ 



2 Notice the difference between Fig. 14 in Paper I and Fig.|2]in the present 
paper. All plots in Paper I and most of the ones in this paper present data 
averaged over 5 orbital periods, however in the case of model M2 the torque 



© 0000 RAS, MNRAS 000,[T][l8] 



4 A. Peplinski et al. 




3.1 Orbital evolution 

Figure [T] shows the evolution of several key quantities for models 
Ml and M2 (curve 1 and 2 respectively). The upper left panel dis- 
plays the evolution of the planet's semi-major axis. Both models 
show two phases of migration. The first phase is one of rapid, out- 
ward directed migration, lasting for about 30.7 and 34.4 orbits in 
Ml and M2 respectively. It is followed by a phase of inward mi- 
gration, which is slower for Ml than for M2. In Ml the effective 
planet mass grows rapidly during the first phase of migration, and 
the planet, after reaching a w 1.72, becomes locked in the disc in a 
type II like migration. In M2 this growth is artificially inhibited, and 
the planet migrates further, reaching a maximum semi-major axis 
at a « 2.26. At this position the planet starts an inward-directed 
rapid migration phase. 

In our investigation we concentrate on the first, rapid outward 
directed migration phase and on the reversal of the migration di- 
rection. We will not describe in detail the second, inward migration 
phase, which corresponds to the transition between type III and 
type II migration regime in Ml, and to the rapid inward migration 
in M2. 

The time evolution of the effective planet mass A/p and the 
mass of the gas inside the Hill sphere A/h are shown in Fig. [2] (see 
also Paper I, Fig. 14, but notice the change of the labels). In both 
models the planet accumulates gas during the outward migration 
phase. In Ml A/h grows faster than in M2, and reaches HM^ at 
30 orbits. This is a considerable fraction of the original mass of the 
protoplanetary disc, which contained 29A/i ) _ within r — 2. In the 
slow inward migration phase the planet's mass is almost constant 
and A/p « IOMl. . Comparing A/p and Mh we can see that most 
of the mass is actually accumulated inside the planet's gravitational 
smoothing length. In M2 a smaller amount of gas enters the Hill 
sphere and a maximum of Mh ~ 5Mi^_ is reached at t = 40 orbits. 

This growth of A/p is possible due to the use of a locally 
isothermal equation of state with a constant circumplanetary disc 
aspect ratio hp. We study the effects of varying h p below in 
Sect. 14.21 Note that in Ml there is no mass outflow from the Roche 
lobe during the slow migration phase, such as was seen in the in- 
ward migration case described in Paper II, Sect. 4.3. In this respect, 

and the migration rate oscillate strongly in time and we instead averaged 
over 10 orbits. 



Ml is more similar to the inward migration simulation with a lower 
h p = 0.3 (see Paper I, Sect. 5.1.2). 

The evolution of the migration rate d is presented in mid- 
dle left panel of Fig. Q] Initially d grows with the planet's semi- 
major axis for both models and has a maximum at f ~ 20 orbits 
(d « 0.0062) and t » 23.3 orbits (d » 0.011) for Ml and M2 
respectively. In Ml d has a lower value due to the increase of the 
planet's inertia. Later d drops in a similar way for both models, 
however in M2 the planet goes through a short period of strong 
oscillations of the migration rate around t « 28 orbits. The physi- 
cal mechanisms behind the drop in migration rates are addressed in 
Sects. EHl and|5l 

During the slow inward migration phase of model Mid oscil- 
lates around —3 x 10~ 4 . The oscillations are slightly more visible 
after about 60 orbits after the planet's eccentricity e has increased. 
As we saw before, in M2 the planet enters a rapid inward migration 
phase. 

The plot of the non-dimensional migration rate Z presented 
in middle right panel helps us understand the difference between 
the two simulations. As explained in Paper H, the non-dimensional 
migration rate is defined as a ratio of the planet's migration rate 
d and the so-called fast migration speed df, which is given by the 
ratio of the half width of the horseshoe region x s and the libration 
timescale Tiib dArtvmowicz & Pe plinski 2008). It is equivalent to 
Z — Tiib/T m ig r , where T m i gr = x s /a is the migration time-scale. 
Using Z we can divide type III migration into a fast \Z\ > 1 and 
a slow \Z\ < 1 migration regime. Note that this division is not 
identical to that between type III and type II migration. Type II like 
migration corresponds to the phase of nearly constant \Z\ <C 1. 

The initial value of Z is close to 1 for both Ml and M2. 
However, due to the increase in A/p and the fact that df ~ 
o" 0,5 (Mp) 2 ' 3 , Z never exceeds 1 for model Ml, and the planet 
migrates in the slow migration regime during the entire simula- 
tion. During the outward migration phase Z oscillates around 0.7, 
although both a and d grow during the first 20 orbits. Later Z de- 
creases reaching its final value in the slow inward migration phase. 
For M2, \Z\ > 1 during essentially the entire simulation, so in the 
outward as well as in the inward migration phase the planet mi- 
grates in the fast migration regime. 

As shown in Paper II and below in Sect. 13.21 the value of \Z\ 
correlates with the asymmetry of the co-rotation region, and gap 
opening occurs for \Z\ <C 1. Model Ml thus experiences gap open- 
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ing, and a smooth transition to inward migration, whereas model 
M2 maintains a strong asymmetry in the co-orbital region, without 
any gap opening. 

The upper right panel in Fig. [T] shows the changes of the ec- 
centricity e. In Ml e is approximately constant during the initial 
outward migration phase and does not exceed 1%. At t ~ 53 or- 
bits, the eccentricity starts to rapidly increase. This coincides with 
an increase of the eccentricity of the outer circumstellar disc (outer 
with respect to the planet's orbit). At about 48 orbits the outer disc 
becomes eccentric (lower left panel in Fig.[5} and the torque T starts 
to oscillate (lower left panel in Fig.[T]l. Note that the oscillations of 
r are driven by the gas placed outside the Hill sphere. The oscilla- 
tions of the torque generated by the gas from within the Hill sphere 
are order of magnitude smaller and start to be visible after about 60 
orbits. In the later stages e oscillates around 4%. In M2 the eccen- 
tricity grows during the outward migration phase reaching 3% and 
decreases in the inward migration phase. 

The lower left panel in Fig. Q] presents the time evolution of 
the total torque exerted by the gas on the planet. T grows during the 
rapid outward migration phase in both models, however in Ml there 
is a short (~ 10 orbits) period of decreasing torque. The torque has 
a maximum at about 20 orbits and its value, unlike the maximum 
value of the migration rate, is similar for both models. Later on F 
decreases rapidly and the planets enter the inward migration phase. 
There is visible a short period of strong oscillations of F at about 
28 orbits in M2. As noted above, F starts to oscillate in Ml at 48 
orbits, when the eccentricity of the outer disc starts to grow. M2 
also experiences a short period of torque oscillations around t ~ 28 
orbits. 

In Paper II Sect. 4.1 we described the relation between the 
averaged torque and the non-dimensional migration rate Z for the 
inward migration case. |T| was shown to grow linearly with \Z\ in 
the slow migration regime, and to decrease slowly in the fast migra- 
tion regime. The trend for the outward directed migration is similar, 
however there are also some differences. As we saw above, models 
Ml and M2 show very different values for Z, and it we can study 
the relation between V and Z in slow (Z < 1) migration regime by 
looking at Ml, and in the fast (Z > 1) migration regime using M2. 
We show this relation in the lower right panel in Fig.Q] In the slow 
migration regime (Ml) the torque evolution is strongly influenced 
by the mass accumulation in the circumplanetary disc and shows 
two different stages. In the first one (starting at Z ~ 1) V grows 
reaching the maximum at Z ~ 0.6. The second stage with V ~ Z 
corresponds to a slowing down of migration and agrees with the 
results for the inward directed migration. Model M2 (always in the 
fast migration regime) shows two stages too. The first one qualita- 
tively agrees with the results for the inward directed migration with 
F oscillating between 0.003 and 0.0045, and being less dependent 
on Z. In the second stage the planet reverts rapidly its direction of 
migration and F decays quickly with Z. 

3.2 Flow structure in the co-orbital region 

In Paper II we discussed the relation between the time-averaged 
torque and the asymmetry of the horse shoe region for the inward 
migration case. We have found the shape of the horseshoe region 
to be dependent on the relation between the migration time scale 
T m i gr and the libration timescale Tub, as expressed by the non- 
dimensional migration rate Z. In the slow migration regime the 
horse shoe region extends almost the entire 27r angle of the co- 
rotation region, and is fairly symmetric. The co-orbital flow is lim- 
ited to a small region near to the planet. In the fast migration regime 



the horseshoe region shrinks to a single tadpole-like region of lim- 
ited azimuthal extent, disconnecting from the planet on one side 
and the co-orbital flow takes place over a wide azimuthal range. 

A similar dependency is observed for the outward migration 
case. However, the fact that the volume (and mass) of the co-orbital 
region grows with a during outward migration produces important 
differences. The relation between the shape of the horseshoe region 
and the non-dimensional migration rate for the outward migration 
in the slow and the fast migration regimes is illustrated by models 
Ml and M2, respectively. The results are presented in Figs.|3]and|4] 
The plots show the mass fraction of the gas that initially was placed 
in the corotational region ai n it — 2Ru < r s < a + 2Rh, where 
r s is here the distance to the star and ainit is the initial planet's 
semi-major axis. The mass fraction goes from (dark blue; region 
occupied by the gas from the inner or outer disc only) up to 1 (red; 
region occupied by the gas from corotation only). 

We first consider the slow migration regime (Ml, Fig. [3). In 
this case, unlike the inward directed migration, most of the gas 
placed initially in the co-orbital region is captured by the planet 
in the tadpole-shaped horseshoe region around the generalised li- 
bration point G^j and migrates together with the planet. The top 
left and a top right panels present the outward migration phase for 
Z « 0.6 (t = 20 orbits) and Z » 0.1 (t = 30 orbits) respec- 
tively. In the first panel the horseshoe region is visibly asymmet- 
ric, with only one tadpole being present. The azimuthal extent of 
the horseshoe region is about 1.77T and the co-orbital flow is a rel- 
atively broad stream, giving a large positive value for the torque 
r. In the second panel both tadpoles are present and the asym- 
metry of the horseshoe region is starting to disappear. As the az- 
imuthal extent of the horse shoe region grows, that of the co-orbital 
flow shrinks. This reduces the co-orbital torque Ten (and the total 
torque F which is dominated by the co-orbital torque) consider- 
ably. The co-orbital flow keeps decreasing with decreasing Z and 
disappears in the subsequent inward slow migration phase, when 
the horseshoe region covers the entire azimuthal range (lower left 
panel). For \Z\ <C 1 the planet's radial motion is slow and a gap 
has started to be cleared as the planet makes the transition to type 
II migration. The last panel displays the mass fraction at t — 70 
after the planet has almost opened a gap in the strongly eccentric 
circumstellar disc. 

During the outward migration phase with Z ~ 1 (upper left 
panel) the gas captured in the horseshoe region does not mix with 
the co-orbital flow. The situation changes for the slower migration. 
The horseshoe region increases its volume with decreasing Z and 
the planet captures part of the gas from the co-orbital flow starting 
a strong mixing of the gas. 

This behaviour of model Ml is similar to the scenario de- 
scribed in Paper II, where we also found a smooth transition be- 
tween the fast and slow migration regimes corresponding to an in- 
crease in the extent and symmetry of the horse shoe region, and a 
shrinking of the co-orbital flow. Note however that the reasons are 
different: decrease of the co-orbital mass for the inward migration 
case and increase of Mp for Ml. 

Model M2 shows a different evolution, see Fig. [4] As in Ml 
most of the gas initially placed in the co-orbital region is cap- 
tured by the planet in the tadpole- shaped horseshoe region around 



3 In Paper II we introduced the concept of the generalised Lagrangian 
points G4 and G5, which for a slowly migrating planet correspond to L4 
and L5 libration points. This distinction is necessary, since for a rapidly 
migrating planet these libration points move in a direction of a planet. 
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Figure 3. The mass fraction of the gas that initially was placed in the corotational region for model Ml . The mass fraction goes from (region occupied by the 
gas from the inner or outer disc only) up to 1 (region occupied by the gas from corotation only). The top left (t = 20 orbits Z ss 0.6) and a top right (t = 30 
orbits Z ps 0.1) panels present the outward migration phase. The lower left and lower right panels show the slow migration phase \Z\ <C 1 at t = 50 and 70 
orbits respectively. The upper row and the lower left panel present the disc for the planet's eccentricity e < 0.01. After about 50 orbits both the planet's and 
the outer disc eccentricity rapidly grow giving a strongly asymmetric flow. It is presented in the lower right panel. All plots are done in co-moving reference 
frame with the planet at (x, y) equal (1.42, 0), (1.71, 0), (1.66, 0) and (1.62, 0) respectively. 



the G5 libration point and migrates together with the planet. The 
top left and a top right panels present the migrating planet dur- 
ing the phase of increasing Z, with Z ~ 1.7 (i = 10 orbits) 
and Z ~ 4.5 (t — 25.3 orbits). At both times the horseshoe re- 
gion is a single tadpole-like region, and the co-orbital flow has a 
wide azimuthal range. The azimuthal extent of the horseshoe region 
clearly decreases with Z and is about 1.7-7T and 0.5-7T for Z = 1.7 
and Z = 4.5 respectively. On the other hand the radial extent of 
this region increases, since x s ~ a, and the amount of mass in 
the co-orbital region grows during the migration making outward 
migration potentially a self-accelerating, run-away process. In Ml 
the run-away character is suppressed due to the rapid grow of the 
planet's inertia. But even in M2, where the planet mass is constant, 
we do not find run-away migration. The reason is that for the given 
disc structure there is a critical value of Z, for which the horse- 
shoe region disconnects from the planet, leading to a decrease of 
the co-orbital torque. This process is shown in the lower left and 
lower right panels. Here the planet is slowing down its migration 
at i = 31.6 (2 RJ 2) and 36 (Z » 0) orbits respectively. The 



pink lines are the flow lines at the border of the horseshoe region 
and show that the gas marked by the red colour does not fill the 
whole horseshoe region. After the initial horseshoe region discon- 
nects from the planet causing a drop in For, the azimuthal extent 
of the horseshoe region grows rapidly as the migration rate drops. 
As a result a large amount of gas from the co-orbital flow is cap- 
tured by the planet in the new horseshoe region, stopping the mass 
transfer between outer and inner disc which was responsible for the 
rapid migration. We find the stream-lines in the co-orbital region to 
change rapidly and no smooth transition from outward to inward 
migration occurs. This behaviour is similar to the stopping of in- 
ward migration at a disc edge (Paper II, Sect. 6). We will discuss 
this in more detail in Sect. [5] However we will encounter these two 
stopping mechanisms below in other simulations. We will refer to 
them as stopping mechanism Ml (for stopping due to mass accu- 
mulation), and M2 (for loss of the horseshoe region). 

Model M2 does not show not much mixing between the gas 
originally captured in the horseshoe region and the co-orbital flow. 
Note that during the entire phase of fast migration a small amount 



© 0000 RAS, MNRAS OOO.mfTsI 



Type III migration; Outward migration 7 



simulation tirne= 6.2&33E+G1 Jyrt. 



simulation tirne= I Jyrt. 





0.2 0.4 0,« 0,« I 

2 

Figure 4. The mass fraction of the gas that initially was placed in the corotational region for model M2. The top left (t = 10 orbits Z 1.7) and a top right 
(t = 25.3 orbits Z PS 4.5) panels present the migrating planet increasing the non-dimensional migration rate. The lower left and lower right panels show the 
planet in the stage of slowing down the migration at t = 31.6 (Z fti 2) and 36 (Z ?s 0) orbits respectively. The pink lines show the flow lines at the border of 
the horseshoe region. All plots are done in co-moving reference frame with the planet at (x, y) equal (1.1, 0), (1.87, 0), (2.1, 0) and (2.25, 0) respectively. 



of the gas is leaving the horseshoe region, flowing along the bow 
shock in the inner disc. This is visible in the plots as a light blue 
spiral. 



3.3 Surface density 

The global surface density distribution is presented in Fig. [5] The 
left column shows the surface density profiles for model Ml at 
t = 20 (upper), 30 (middle) and 70 (lower) orbits. These corre- 
spond to upper left, upper right and lower right panels in Fig. [3] 
The first plot shows the planet migrating outward with Z « 0.6. 
The shape of the low density region in the co-orbital region agrees 
mostly with the shape of the horseshoe region, but fills the whole 
azimuthal range and the first signs of gap creation are visible. 
This is more prominent at t — 30 orbits (Z « 0.1), when the 
horseshoe region fills the whole corotation region and captures 
gas from the co-orbital flow. There is a small asymmetry between 
the Lagrangian points L4 and L5, which does not grow until the 
planet's eccentricity increases rapidly. After about 50 orbits both 
the planet's and the outer disc eccentricity quickly grow and the 



interaction of the disc with the outer boundary starts to be impor- 
tant. After 120 orbits about 15% (14Mq, ) of the disc mass leaves 
the computational domain through the outer boundary. The surface 
density for the eccentric outer disc is presented in the last plot. 

The right column presents the surface density for model M2 
at times 25.3 (upper), 31.6 (middle) and 36 (lower) orbits, corre- 
sponding to the upper left, upper right and lower right panels in 
Fig. [4] There is a visible asymmetry in the co-orbital region in all 
the plots, but only in the first and second plot is this related to the 
co-orbital flow. These two plots show the planet migrating in the 
fast migration regime with Z ~ 4.5 (first plot) and Z » 2 (second 
plot). The first one corresponds to the maximum value of Z (and 
d) and shows the smallest extent of the horseshoe region during the 
entire simulation. The low density region represents the gas that 
was initially captured by the planet and covers the whole horseshoe 
region. The strong co-orbital flow is visible. The second plot rep- 
resents the phase when the planet has reduced its migration rate. 
Here the horseshoe region has grown and captured gas from the 
co-orbital flow. The low density region does not cover the whole 
horseshoe region and only represents the position of the gas that 
was initially captured by the planet. It has sharp edges as no mixing 
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Figure 5. Global surface density for the outward migration case. Left and right panels correspond to models Ml and M2 respectively. The upper left and 
middle left panels present the outward migration phase at t = 20 (Z 0.6)and 30 (Z ta 0.1) orbits and the planet's eccentricity e < 0.01. The lower left 
panel shows the inward, slow migration phase Z\ <C 1 at 70 orbits and e Si 0.04. The top right panel presents the stage of migration for the maximum value 
of the non-dimensional migration rate (t = 25.3 orbits Z m 4.5). The middle right and lower right panels show the planet in the stage of slowing down the 
migration at t = 31.6 (Z Si 2) and 36 (Z ss 0) orbits respectively. All plots are done in co-moving reference frame with the planet. The colour scale is 
logarithmic and differs between Ml and M2. 
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Figure 6. The surface density profiles for models Ml (left column) and M2 (right column). The upper row shows the azimuthal average of the surface density 
E. Upper left panel presents it at time t equal 0, 10, 20, 30 and 40 orbits (curves 1, 2, 3, 4 and 5; Z equal 0, 0.7, 0.6, 0.1 and —0.02 respectively ). Upper 
right gives the density profile at 0, 10, 20, 25.2 and 36 orbits (curves 1, 2, 3, 4 and 5; Z equal 0, 1.7, 2.7, 4.5 and respectively). The planet position is visible 
as a strong spike. The lower row presents the initial density profile (curve 1 ), the azimuthal average of the surface density (curve 2) and the surface density 
cuts through the libration points L5 (curve 3) an L4 (curve 4). Lower left and lower right panels correspond to time equal 20 (Z m 0.6) and 25.2 (Z as 4.5) 
orbits respectively. 



has happened. The last plot shows the planet reversing its direction 
of migration with Z ~ 0. There is a visible asymmetry in the coro- 
tation, even though the co-orbital flow is very weak here. This is 
caused by the fact that the transition between the outward and in- 
ward migration takes place rapidly, and the gas does not have time 
to adjust to the new configuration. For the same reason the planet 
does not open a gap, even though it temporarily migrates in the 
slow migration regime. At this stage the low density region starts 
to disappear, and it is completely gone during the inward migration 
phase. 

These processes are further illustrated in the upper row in 
Fig. [6] The upper left panel presents the azimuthal average of the 
surface density E over 2tt for model Ml at times t = 0, 10, 20, 30 
and 40 orbits (curves 1, 2, 3, 4 and 5). The first curve shows the ini- 
tial conditions, and curves 2, 3 and 4 present the outward migration 
phase. Curve 5 corresponds to the inward migrating planet inside a 
gap. The planet in fact starts to clear a gap shortly after the migra- 
tion starts, however the co-orbital flow is strong enough to fill the 
inner region of the disc, increasing the density there considerably. 
The upper right panel presents the azimuthal average of the surface 
density for model M2. Curve 1 shows the initial conditions, and the 
curves 2, 3 and 4 correspond to the outward migrating planet. Curve 
5 shows the average density profile for the planet at the transition 



between outward and inward directed migration Z m 0. During the 
whole simulation the planet is unable to open a gap, even during 
the time when Z ~ 0. Like in model Ml the co-orbital flow is 
strong enough to fill the inner disc and the final disc density profile 
is relatively constant. 

The lowest row in Fig. [6] illustrates the density asymmetry in 
the co-orbital region. It shows the initial density profile (curve 1), 
the azimuthal average of the surface density (curve 2) and the sur- 
face density cuts through the libration points L5 (curve 3) an L4 
(curve 4). The cut through L5 gives the approximate width of the 
co-orbital region and the density E g inside the horseshoe region. 
Similarly the cut through L4 gives the approximate value of the 
density of the gas crossing the co-orbital region E s . As explained 
in Paper II, the difference between E s and E g is the main reason 
for type III migration. 

The lower left panel presents Ml at 20 orbits during the out- 
ward migration phase {Z ~ 0.6). Since the planet is in the slow mi- 
gration regime and gap opening has already started, the asymmetry 
between L4 and L5 is limited to the outer edge of the co-orbital 
region. In contrast, model M2 (lower right panel) has the planet 
migrating in the fast migration regime and the asymmetry covers 
the whole co-orbital region. In this case the planet carries along the 
low density horseshoe region with well defined sharp edges and an 
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Figure 7. Results of the simulations for the different planet's masses. Curves 1, 2 and 3 correspond to the initial planet mass Mp equal 0.0007, 0.001 (standard 
case Ml) and 0.0013 respectively. Upper left and upper right panels show the planet's semi-major axis a and the migration rate a during the orbital evolution. 
The non-dimensional migration rate Z as a function of the planet's position in the disc and the gravitational planet mass Mp are presented on the lower left 
and lower right panels respectively. 



almost static value of the surface density E g « 0.0005. The value 
for E s is close to the initial density profile at the current planet's 
position. 



4 DEPENDENCE ON THE SIMULATION PARAMETERS 

After having described the two standard cases of outward migra- 
tion, we now discuss the dependence of the migration process on 
various simulation parameters. 

4.1 Planet's mass 

The first important parameter is the planet mass. In our investi- 
gation we concentrate on the orbital evolution of giant planets. In 
Papers I and II we showed that in the fast migration regime, planet 
mass does not influence the migration significantly, whereas in the 
slow migration regime the migration is strongly dependent on Mp . 
The relation between Mp and the planet's migration for outward 
migration is presented in Fig. [7] We ran models with initial masses 
M P equal 0.0007, 0.001 (standard case Ml) and 0.0013 (curves 
1, 2 and 3 respectively). The effective planet mass was allowed to 
increase, Mp = Mp. 

The upper left and upper right panels show the planet's semi- 
major axis a and the migration rate a. There is a clear difference 



between the evolution of the two first models and the last one, in- 
dicating a sharp border in parameter space for which the outward 
migration is allowed. In the last model the amount of mass in the 
corotation is too small to support the initial impulse for outward mi- 
gration and the planet transits to an inward migration case at a ~ 1. 
Note that even though the initial migration rate a is similar for all 
the models, the non-dimensional migration rate (lower left panel) 
differs since Z ~ M^ 2 ^ 3 . We find that to start outward migration 
the initial average Z should be bigger or close to 1, and the last 
model does not satisfy this criterion having an (extrapolated) initial 
value of Z ~ 0.9. Mass accumulation does not play a role here, as 
the increase of Mp is relatively modest (~ 2.5M], after 10 orbits) 
and does not influence the migration in a major way. 

The two lower mass models satisfy the criterion stated above, 
with extrapolated values of Z ~ 1.15 and Z ~ 1.05, respectively. 
During the first 25 orbits both systems evolve almost exactly the 
same way, even though the migration takes place in the slow mi- 
gration regime. This can be understood from the fact that the non- 
dimensional migration rate is relatively close to 1 (Z > 0.65) and 
the migration is mostly determined by the increase of the gravita- 
tional planet's mass, which grows in a similar way in both models 
(reaching about lOM^ ). Unlike the inward migration case (Paper II 
Sect. 5.1), the amount of mass accumulated by these planets is al- 
most independent of its initial mass. The difference is relatively 
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small (of order of Q.bMi.) and does not play the role while a is 
increasing. 

The evolution starts to differ after Z drops below 0.6 and a 
starts to decrease. The migration rate drops faster for the model 
with the larger initial planet mass and the planet moves from the 
outward to inward migration at about 30 orbits and a ~ 1.7. For 
the lighter planet this transition takes place at about 35 orbits and 
a « 1.9. 

The relation between planet mass and orbital evolution is vis- 
ible on the lower left panel, where the non-dimensional migration 
rate Z is plotted as a function of the planet's position in the disc. 
The increase of planet mass, until Mp reaches the critical mass 
for which outward migration is no longer possible, causes Z to di- 
minish smoothly and moves the curve to the right. In this case the 
region of the disc where the rapid migration can take place shrinks 
with increasing Mp and the planet stops closer to the star. 



4.2 Effects of circumplanetary disc aspect ratio 

The next parameter we investigate is the circumplanetary disc as- 
pect ratio h p . We performed a set of three simulations with dif- 
ferent values of h p and Mp = Mp. The results are presented in 
Fig. [8] We use the same set of plots as in Fig. [7] the evolution of 
the planet's semi-major axis a and the migration rate a in the upper 
panels, and the non-dimensional migration rate Z as a function of 
the planet's position in the disc together with the evolution of the 



effective planet mass in the lower panels. Curves 1, 2 and 3 corre- 
spond to hp = 0.4 (standard case Ml), 0.5 and 0.6 respectively. 

The planet's mass Mp is strongly dependent on the circum- 
planetary disc aspect ratio and decreases with increasing h p . In the 
first model with h p = 0.4, the gravitational planet's mass grows 
up to lOM-i,. causing the planet to migrate in the slow migration 
regime. In the second model with h p = 0.5, Mp grows during the 
whole outward migration stage reaching 2A/-2 f , however it stays 
very low during the first 10 orbits allowing the planet to migrate in 
the fast migration regime. A similar behaviour is seen for h p — 0.6, 
but in this case Mp is almost constant. 

There is an important difference between the first (h p = 0.4) 
and the second {h p = 0.5) model. In the first model the outward 
migration stops due due to the rapid increase of the planet's iner- 
tia, but in the second one Mp does not grow enough to reduce Z 
considerably and a increases significantly. This gives a rapid in- 
crease of the non-dimensional migration rate up to Z ~ 3 at which 
point changes in the flow structure in the co-orbital region similar 
to those found for standard case M2, slow down the migration. The 
increase of the circumplanetary disc aspect ratio allows the planet 
to migrate faster and to larger radii. The planet reaches a ~ 1.72 
at t ~ 21 orbits and a ~ 1.9 at t ~ 24 orbits in the first and the 
second model respectively. 

A further increase of h p does not influence the migration sig- 
nificantly, and the evolution of the second and the third model is 
similar. However, increasing h p implies lowering Mp, giving a 
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Figure 9. Results of the simulations for the different total disc masses for Mp = Mp . Curves 1 , 2 and 3 correspond to the disc to the primary mass ratio /ip> 
equal 0.0075, 0.01 (standard case Ml) and 0.0125 respectively. The panels are the same as in Fig. [7] 



faster increase of Z. That is why in the model with h p — 0.6 the 
flow inside the Roche lobe changes its structure for a lower value 
of a and the planet changes to inward migration at t ~ 23 orbits 
reaching a maximum orbit a ~ 1.7. In the last model the migration 
rate grows up to 0.09, but in the second model a reaches 0.11. 

Summarising, as in the inward migration case the temperature 
profile influences the planet's orbital evolution mostly by limiting 
the amount of mass in the Hill sphere. h p ^ 0.4 gives a rapid 
grow of Afp and the planet stops due to the increase of its inertia. 
Higher values for h p reduce the effective planet mass and lead to 
the other stopping mechanism (M2), where the flow in the planet's 
vicinity changes structure for some critical value of Z. In this case 
the lower gravitational planet mass (corresponding to higher h p ) 
allows to reach this critical value of Z faster and stop outward mi- 
gration closer to the star. 

The reader has to keep in mind that we present the results 
of two-dimensional simulations with h p constant in time, an im- 
portant simplification. In the case of strong mass inflow into the 
planet's vicinity, h p should be a function of the gas density in the 
circumplanetary disc, and the possible vertical motions of the gas 
should be taken into account. 



4.3 Dependence on the total disc mass and the initial density 
profile 

Above we discussed the effects of parameters connected to the 
planet and its circumplanetary disc. In this section we consider 



parameters describing the global circumstellar disc: the total disc 
mass given by the disc to the primary mass ratio /id , and the initial 
density profile described by the exponent qe- 



4.3.1 Total disc mass 

To describe the effect of changing the total disc mass we performed 
two sets of simulations: one set for the gravitational planet mass 
increased by the mass within the planet's gravitational softening 
Mp — Mp (similar to Ml) and one set for the constant planet's 
mass Afp = Mp (similar to M2). In both sets Mp = 0.001 and 
a E = -0.5. 

The first set contains three models with no equal 0.0075, 0.01 
(standard case Ml) and 0.0125 (curves 1, 2 and 3 in Fig.[9j(. We see 
that the results are similar to those of the simulations performed 
for different planet masses (Sect 14. ll and Fig. [7}. As in the case 
of inward migration (Paper II Sect. 5.3), an increase of the disc 
mass acts in a similar way as decreasing the planet mass. There are 
however also some differences. 

In the simulation with the lowest disc mass the amount of the 
mass within corotation is not sufficient to support type III migra- 
tion, and the planets stops outward migration after 14 orbits reach- 
ing a fa 0.96. In this case Z « 0.85 at the start of the simulation, 
and the Z > 1 criterion for outward migration is not satisfied. This 
illustrates again the sharp border in parameter space for which out- 
ward type III migration is allowed. The combination of planet and 
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Figure 10. Results of the simulations for the different total disc masses for Afp = Afp = const. Curves 1 and 2 correspond to the disc to the primary mass 
ratio ud equal 0.0075 and 0.01 (standard case M2), respectively. The panels are the same as in Fig. [7] however the lower right panel shows the changes of 
the mass inside the Hill sphere during evolution instead of effective planet's mass. 



disc masses gives a single parameter defining the disc region where 
the type III migration is allowed. 

In the second and the third model the condition for the out- 
ward migration is satisfied with the initial averaged Z equal 1.05 
and 1.1 respectively. The evolution of these models is similar, how- 
ever the migration rate for the most massive disc is slightly higher 
and the planet reaches a ~ 1.73 at 27 orbits. In the second model 
the outward migration stops at 32 orbits and a ~ 1.7. The last 
model shows the fastest increase of the gravitational planet mass 
too, giving Afp ~ lSMq, after 27 orbits. In the second model 
Afp ~ 10Mq^ at 30 orbits. This reduces the difference between 
the migration rate d, and the plot of the non-dimensional migra- 
tion rate as a function of radius (lower right panel) shows a very 
similar evolution for both systems. This is different from the mod- 
els with different planet masses (where the lighter planet can travel 
further away), and the inward migration case (where no modifies 
the planet's migration in a significant way). 

In the second set of simulations (Afp = Afp) we have two 
models with /id equal 0.0075 and 0.01. The results are presented in 
Fig-Co] Curves 1 and 2 correspond to the first and the second model 
respectively. In a more massive disc the planet migrates faster and 
has a higher mass accumulation rate, but the maximum value of 
the mass within the Hill sphere A/h is similar for both models and 
equals about 5M-^ . Similarly the migration rate increases faster for 
the more massive disc, but both models reach the same maximum 
value, a ~ 0.01 at 43 and 25 orbits for model 1 and 2 respectively. 



After the migration rate starts to decrease, A/h is still growing and 
reaches its maximum at 55 and 40 orbits. However, the transition 
between outward and inward directed migration takes place at 50 
(a ~ 2.6) and 35 (a ~ 2.25) orbits. So, although the more massive 
disc allows the planet to travel faster, the outward migration phase 
lasts shorter and the planet actually stops closer to the star. This is 
caused by the nature of the stopping mechanism which here is M2- 
like. In both models the rapid changes start upon reaching Z ~ 4, a 
value only weakly dependent on /id • The more massive disc shows 
a stronger increase of Z and reaches this critical value sooner. 

4.3.2 Initial density profile 

The final parameter we consider is the density exponent as - Since 
the evolution of the models with Afp — Afp is dominated by the 
mass accumulation and only weakly depends on fj,u , we performed 
a series of simulations with constant Afp = Afp, /id = 0.0075 
and «s equal 0.0, —0.5, —1 and —1.5 (curves 1, 2, 3 and 4 in 
Fig. lilt . The density gradient at the initial planet position is the 
same in all presented models, giving the same initial value of the 
migration rate a. 

Once again this set of simulations shows the sharp boundary in 
parameter space for which outward migration is allowed. Although 
the initial condition for outward migration Z > 1 is satisfied for 
all the models, in the last simulation (as = —1.5, curve 4) the 
mass within corotation is too small to support the migration and the 
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Figure 11. Results of the simulations for the different initial disc profiles with Afp 
— 1 and —1.5 respectively. The panels are the same as in Fig. 1101 
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Mp = 0.001. Curves 1, 2, 3 and 4 correspond to equal 0.0, —0.5, 



planet switches to inward migration after about 10 orbits reaching 
a ~ 1. The mass within its Hill sphere does not exceed 0.7M%_ . 

The other models migrate in the fast migration regime, how- 
ever the migration rate depends strongly on as. The smallest gra- 
dient in the density profile corresponds to the biggest mass in the 
corotation region and allows for the fastest increase of d, reach- 
ing 0.015 after 30 orbits. This gives a maximum value of Z ss 5.5. 
Similarly this model shows the quickest increase of the mass within 
the Hill sphere giving Mh ft! 5Mi, at about 40 orbits. The transi- 
tion between outward and inward migration takes place at about 37 
orbits and a w 2.4. 

Lower as means lower mass within the corotation region, and 
thus slower migration. That is why the migration rate grows slower 
in the second model reaching 0.01 after about 46 orbits, and re- 
mains approximately constant with a « 0.005 between 20 and 50 
orbits in the third model. Although a is almost constant in model 
3, Z grows in both models reaching about 4 and 2.5 in the second 
and the third model respectively. The lower value of Z allows the 
planets to travel outward for a longer time, and only switch to the 
inward migration phase after about 51 (a « 2.6) and 56 (o « 2.3) 
in model 2 and 3 respectively. The rate of the mass accumulation 
in the Hill sphere decreases with decreasing as, however the final 
mass in the first and the second model is quite similar. In the third 
model Mh reaches almost 4Mi_. 

One thing that this series of simulation shows is that the criti- 
cal value of Z for stopping mechanism M2, is not universal, but ac- 



tually depends on the exponent in the density profile and decreases 
with decreasing as- 



5 STOPPING TYPE III MIGRATION 

In Paper II we described two possible mechanisms of stopping 
type III migration (or more precisely, switching it to type II migra- 
tion). Inward migration can be slowed down by the simple geomet- 
rical effect of the shrinking of the co-orbital region with decreasing 
semi-major axis, or by the interaction with a steep density gradi- 
ent, such as a disc edge. The second mechanism acts in the same 
way for outward directed migration, but the first one would make 
outward type III migration a self-accelerating process. In this case 
an additional mechanism is necessary to limit the non-dimensional 
migration rate We found two possible mechanisms to achieve 
this. 

The first one is the rapid increase of the planet's inertia, which 
lowers Z and forces the migration to proceed in the slow migration 
regime (Z < 1). The example is our standard case Ml, in which 
the effective planet mass grows from 1 up to about 10Ma during 
about 30 orbits (~ 350 years). Such a rapid mass accumulation 
seems unrealistic and could be expected to lead to an increase of 



4 In most of the simulations a grows during the outward migration phase 
even though Z decreases; see middle row in Fig.[TJ 
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Figure 12. The non-averaged torque exerted by the gas on the planet (upper left panel) and the details of the gas flow in the Roche lobe for model with 
hp = 0.6. Upper right, lower left and right panels show the surface density and the flow lines in the planet's vicinity after 12.2, 13 and 13.2 orbits. The colour 
scale is logarithmic. The plotted domain is square region of the size of 6-Rh ■ The flow lines close to the border of different regions are made more visible. 



the disc aspect ratio h p , which would reduce the mass accumula- 
tion rate. Our disc model does not allow for time-dependent values 
of hp, but as the growth of the rapidly migrating planet is a highly 
dynamical process in which a large amount of mass may enter the 
Roche lobe and interact with the gaseous planet's envelope, a much 
more sophisticated method would be needed to model this realisti- 
cally. 

Note that even though such a rapid increase of the effective 
planet's mass is more likely to happen in outward migration, we 
even found it in some cases of inward migration (see the model 
with hp — 0.3 described in Paper I Sect. 5.1.2). 

The second stopping mechanism is encountered in models in 
which the growth of the effective planet mass remains limitetjfl 
and the example is our standard case M2. Here both a and Z grow 
quickly during the migration until a critical value of Z is reached 
and the flow structure in the co-orbital region changes dramati- 
cally. The reason for such a rapid change are strong oscillations 
of the flow in the planet's vicinity. They are clearly visible in the 

5 In the model with hp = 0.5 the final effective planet's mass is of the 
order of 2M % ; see Fig.[8] 



simulation with hp — 0.6 (see Fig. [8). In this model the outward 
directed rapid migration lasts for about 23 orbits and the planet 
reaches a ~ 1.7. This simulation has the important advantage that 
the relatively high value of the circumplanetary disc aspect ratio 
suppresses the highest frequency oscillations, and consequently it 
is easier to study the non-time-averaged torque V. 

The results are presented in Fig. [12] The upper left plot shows 
the time evolution of the non-averaged torque. There are oscilla- 
tions of r that start at about 10 orbits (Z k, 2) and grow in time. 
These oscillations are caused by cyclic changes of the shape of the 
circumplanetary disc. One of these cycles starting at about 12 or- 
bits, is presented in the upper right and lower panels. In Paper II 
we described the flow structure in the planet's vicinity for differ- 
ent stages of migration, showing that in the fast migration regime 
the circumplanetary disc is strongly asymmetric (Paper II, upper 
left panel in Fig. 4) and gives a significant contribution to the total 
torque (Paper II, Fig. 5). A similar situation can be seen in the up- 
per right panel in Fig. [12] corresponding to a value of Y close to the 
maximum of the cycle. The circumplanetary disc has a triangular 
shape compressed at the side of the co-orbital flow and stretched 
at the side of the horseshoe region. In a later phase, the orbits in 
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Figure 13. The eccentricity e as a function of the migration rate a. Left and right panels present the results of the inward and outward directed migration 
respectively. All curves present models with h p = 0.4 and different planet's mass. Curve 1 gives the model with Mp* = Mp = 0.001. Curves 2, 3 and 4 
show the models with Mp* = Mp, and M P equal 0.0007, 0.001 and 0.0013 respectively. 



the circumplanetary disc have stretched even more at the side of 
the bow shock in the outer disc (lower left panel) giving a small in- 
crease of F. In the final phase, some mass from the planet's vicinity 
leaves the Roche lobe flowing along the outer bow shock (lower 
right panel), reducing T significantly. Using the tracer fluid (mark- 
ing the position of gas placed in the initial corotation region) we can 
see that this process involves a big portion of the gas from the co- 
orbital flow entering the Roche lobe and being redirected back to 
the outer disc and not accumulating in the planet's proximity. This 
strongly perturbs the outer bow shock and giving the oscillations in 

r. 

We find that these oscillations grow in time, until they become 
so strong that they destroy the regular co-orbital flow, starting a 
period of strong oscillations in T and a. These strong oscillations of 
the migration rate open the flow-lines in the co-orbital region and a 
large amount of gas from the former co-orbital flow is captured into 
the new horseshoe region. Such a short phase of strong oscillations 
of averaged T and a is also visible in the standard case M2 (Fig.[T](. 

We note that we found a similar process (the corrugation of 
the inner bow shock by the gas from the co-orbital flow being redi- 
rected back to the inner disc) in models for inward migration. How- 
ever there the oscillations do not appear to grow in time. This dif- 
ference in behaviour may be due to the inherently self-decelerating 
and self-accelerating character of inward and outward migration. 

One may wonder whether this rearrangement of the co-orbital 
region is not caused by numerical effects. Although we cannot fully 
rule out that our disc model influences this process, we still think 
it is a rather robust mechanism. We find it in models with different 
parameters and envelope masses, indicating that it is not our self- 
gr avity correction that plays a r ole here. We also point to the results 
of iMasset & Papaloizoij J2003I) which showed an outward directed 
migration to revert its direction after a relatively short time. 



6 ECCENTRICITY EVOLUTION 

In Paper II we discussed the eccentricity evolution for inward mi- 
grating planets. In this case the eccentricity is damped during the 
rapid migration phase on a time-scale similar to that of migration. 
Moreover we found the damping time-scale to be independent of 
the planet mass. In the case of outward migration we find that the 



eccentricity grows with the migration rate a during the whole out- 
ward migration phase. These results suggest a correlation between 
e and d, as shown in Fig. [T3] This plot shows the results of the 
inward (left panel) and outward (right panel) directed migration 
models for different prescriptions of the effective planet's mass 
(M P * = M P = 0.001 for curve 1 and M P * = M P for curves 
2, 3 and 4) and different initial planet's masses (Mp equal 0.0007, 
0.001 and 0.0013 for curves 2, 3 and 4 respectively). For the out- 
ward directed migration only the rapid migration phase is plotted. 

The eccentricity is approximately constant for |d| < 0.005 
and grows rapidly for bigger absolute values of the migration rate. 
This relation is similar for inward and outward directed migration 
and is weakly dependent on the planet mass. The origin of this re- 
lation is unclear, but it may be caused by the fact that the spiral 
orbit of a rapidly migrating planet cannot be properly represented 
by Keplerian elements. This conclusion is supported by the fact 
that during the rapid migration phase, unlike in the gap formation 
phase, the true anomaly 6 for the migrating planet oscillates in a 
relatively narrow range of angles, but the argument of periastron u) 
ranges from up to 2-7T (Fig.l 14b. 



7 CONCLUSIONS 

In this paper we investigated the outward directed type III migration 
of the high-mass planet embedded in a disc. Using two dimensional 
numerical simulations, performed in a Cartesian coordinate system 
and the inertial reference frame, we studied the orbital evolution of 
planets and the dependency on various numerical parameters. The 
adaptive mesh refinement allowed us to achieve high resolution in- 
side the Roche lobe and study in detail the flow structure in the 
planet's vicinity. To avoid problems with numerical convergence 
we used a modified version of the usual local-isothermal approxi- 
mation, where the temperature depends on the distance to both the 
star and the planet, and also added a correction for the gas self- 
gravity near the planet. 

Both inward and outward type III migration are driven by the 
co-orbital torque. However, the existence of the (inward directed) 
differential Lindblad torque can be expected to break the symmetry 
between the two cases. Indeed, the simulations show that outward 
migration always reverts to inward migration. In fact in some cases 
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Figure 14. Evolution of a true anomaly 9 (left panel) and an argument of periastron ui (right panel) for an inward migrating planet. The rapid migration phase 
lasts here for about 85 orbits and is followed by the gap formation phase. During the rapid migration phase oscillates in a relatively narrow range of angles 
and lu ranges from up to 2n. The situation reverses in the gap formation phase. 



this reversal happens almost instantly, and the initial value of the 
non-dimensional migration rate Z is found to be crucial here. Plac- 
ing the planet at an initial density gradient we force the planet to 
start rapid migration, however the strength of the initial impulse 
depends on the planet mass and the local density profile, and it can 
be measured through the initial averaged value of Z. We found that 
the disc supports outward migration if the initial Z > 1. Since 
Z ~ aAfp _2,/3 the initial Z decreases with the planet mass and 
grows with the disc mass. It means that outward migration is easier 
to start for less massive planets in more massive discs. Note here 
that although initially Z > 1 is required, this does not mean that 
outward type III migration can not later proceed in the slow migra- 
tion regime Z < 1. 

The second important difference between inward and outward 
migration is caused by the relation between the planet's semi-major 
axis a and the volume of the co-orbital region. In Paper II we 
showed that the migration rate can be approximated to be depen- 
dent on the mass of the co-orbital region, and for as > —2 the 
inward rapid migration is a self-decelerating process (the mass of 
the co-orbital region decreases during migration). For the same rea- 
son outward directed migration is a self-accelerating process. The 
migration rate a grows with a during the outward migration phase, 
however the non-dimensional migration rate Z can grow or de- 
crease depending on the effective planet mass. 

We found two different stopping mechanisms for outward type 
III migration. Because of the presence of the Lindblad torques, 
stopping here always implies a reversal of the migration direc- 
tion. The first stopping mechanism occurs when the effective planet 
mass increases enough to make Z fall below 1 so that the migra- 
tion proceeds in the slow migration regime. In this case some as- 
pects of the evolution of the system resemble the evolution of the 
standard case of inward migration described in Paper II. In both 
cases the transition between Z ~ 1 and Z <C 1 is smooth and 
proceed in a similar way, however the reason for slowing down mi- 
gration is different (decrease of the co-orbital mass and increase of 
Mp for the inward and outward migration respectively). The horse- 
shoe region is initially restricted to a single tadpole-like region and 
gradually grows filling the whole co-orbital region and decreasing 
the co-orbital flow. As in the inward migration case the torque is a 
linear function of Z. In this model the migration rate is relatively 



small and the planet starts clearing a gap shortly after the simula- 
tion starts. 

The other mechanism of stopping outward type III migration 
is seen in models where the first mechanism does not operate, and 
the planet keeps migrating in the fast migration regime Z > 1. The 
migration is accelerated due to the growth of the co-orbital mass 
and the asymmetry of the flow in the planet's vicinity grows too, 
leading to corrugations of the external bow shock and oscillations 
in the torque. These oscillations grow in time and finally destroy 
the co-orbital flow, opening the flow lines in the co-orbital region. 
Due to this the planet rapidly stops its outward migration and starts 
migrating inward (also in the fast migration regime). The simula- 
tions show that there is a critical value of Z at which the strong 
oscillations of a start and the planet reverts its migration direction. 
We found this critical value of Z to be weakly dependent on the 
total disc mass fiu, but sensitive to as. For our set of simulations 
this critical value for Z fell in the range 3-6. 

For the second stopping mechanism the whole process of di- 
rection reversion is very brief and the planet does not have time 
to initiate gap clearing. For the same reason the gas from the for- 
mer (before the slowing down of the migration started) horseshoe 
region does not mix with the gas captured in the new horseshoe 
region and can be seen as a low density region in the corotation re- 
gion. Compared to the slow migration regime, the averaged torque 
has a weaker dependency on Z. 

Our simulations show a correlation between eccentricity and 
migration rate, which is almost independent on the planet mass. We 
found the eccentricity to be approximately constant for |d| < 0.005 
and to grow rapidly for bigger absolute values of the migration rate. 
This correlation may be caused by the fact that the spiral orbit of a 
rapidly migrating planet is not well represented by Keplerian ele- 
ments. 

Outward directed type III migration is an extreme process in 
terms of mass accumulation by the planet, and clearly to make 
progress a more sophisticated model for how the planet deals with 
such a massive accretion rate is needed. This would involve both 
improvements in the thermal model used, as well as the inclusion 
of the self-gravity of the gas. However, outward type III migration 
can clearly add interesting aspects to planet formation scenarios, al- 
lowing both rapid growth and a (perhaps only temporary) increase 
of its orbital radius. 
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